Optimal search strategies for hidden targets 
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What is the fastest way of finding a randomly hidden target? This question of general relevance is 
of vital importance for foraging animals. Experimental observations reveal that the search behaviour 
of foragers is generally intermittent: active search phases randomly alternate with phases of fast 
ballistic motion. In this letter, we study the efficiency of this type of two states search strategies, 
by calculating analytically the mean first passage time at the target. We model the perception 
mecanism involved in the active search phase by a diffusive process. In this framework, we show 
that the search strategy is optimal when the average duration of "motion phases" varies like the 
power either 3/5 or 2/3 of the average duration of "search phases", depending on the regime. This 
scaling accounts for experimental data over a wide range of species, which suggests that the kinetics 
of search trajectories is a determining factor optimized by foragers and that the perception activity 
is adequately described by a diffusion process. 

PACS numbers: 



Searching for a randomly located object is one of the 
most frequent tasks of living organisms, be it for ob- 
taining food, a sexual partner or a shelter lj. In these 
examples, the search time is generally a limiting factor 
which has to be optimized for the survival of the species. 
The question of determining the efficiency of a search be- 
haviour is thus a crucial problem of behavioral ecology, 
which has inspired numerous experimental^, H 0,13 
and theoretical^, 0, El works . It is also rel- 
evant to broader domains such as stochastic processes 
theory fill applied mathematics 0] and molecular 
biologv [l4ll5| . 

Anyone who has ever lost his keys knows that instinc- 
tively we adopt an intermittent behaviour combining lo- 
cal scanning phases and relocating phases. Indeed, nu- 
merous studies of foraging behaviour of a broad range of 
animal species show that such an intermittent behaviour 
is commonly observed and that the durations of search 
and displacement phases vary widely Q, 0, 0- The spec- 
trum, which goes from cruise strategy (ex. for large fishes 
that swim continuously such as tuna), to ambush or sit- 
and-wait search, where the forager remains stationary for 
long periods (such as rattlesnake), has never been inter- 
preted quantitatively. The intermittent strategy, often 
referred to as "saltatory" 0, can be understood in- 
tuitively when the targets are "difficult" to detect and 
sparsely distributed, as it is the case for many foragers 
(such as ground foraging birds, lizards, planktivorous 
fish...): since a fast movement is known to significantly 
degrade perception abilities 0,0], the forager must search 
slowly. Then, it has to relocate as fast as possible in or- 
der to explore a previously unscanned space, and search 
slowly again. 

Even though numerous models based on optimiza- 
tion of the net energy gain0 0, 0] predict an optimal 




FIG. 1: A two states model of saltatory search behaviour. 



strategy for foragers, the large number of unknown pa- 
rameters used to model the complexity of the energetic 
constraint, renders a quantitative comparison with ex- 
perimental data difficult. Here, as has already been 
suggested[H 0], we assume that the search time is the 
relevant quantity optimized by the forager in order to ob- 
tain a sufficient daily amount of food and to precede other 
competing foragers. We treat the energy cost only as an 
external constraint that sets the maximal speed of the an- 
imal. We develop a general purely kinetic model of target 
search, which captures the essential features of saltatory 
search behaviour observed for foragers in experiments^, 
when the predator has no information about the prey 
location. 

The central point of our schematic model (see figure 
^1 is that it relies on the explicit description of searching 
trajectories. In particular, as we will show, it permits us 
to elucidate the nature of the search phase. In the fol- 
lowing we assume that the searcher displays alternatively 
two distinct attitudes: 

(1) a search phase, hereafter referred to as phase 1, 
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during which the searcher explores its immediate vicin- 
ity using its sensory organs. As justified below, this lo- 
cal scanning is modelled as a "slow" diffusive movement 
(a continuous random walk with diffusion coefficient D). 
The target is found when this movement reaches the tar- 
get location for the first time. 

(2) a motion phase, referred to as phase 2, during which 
the searcher moves "fast" and is unable to detect a target. 
These repositioning moves are characterized by a ballistic 
motion (at constant velocity v). 

Next, We assume that the searcher randomly switches 
from phase 1 (resp. 2) to phase 2 (resp. 1) with a rate per 
unit time fa (resp. fa), and that the preys are immobile 
and randomly distributed with uniform density. Note 
that the average durations of phases 1 and 2 are then 
l//i and 1/fa. The stochastic nature of the model is 
justified in the case of foragers by the observation of large 
fluctuations of the duration of search and motion phases 
for each animalQ. Indeed, exponential laws are widely 
observed for animal behaviours 0. Our goal is to study 
the optimal search strategy with respect to its kinetic 
mechanism only, i.e. to determine the rates fa, fa which 
minimize the first passage time of the searcher at a target 
location. 

Let us justify briefly our modelling of each phase. It 
is clearly beyond the scope of our model to describe the 
sensory search phase in details, as it involves very com- 
plex biological mechanisms. Our goal here is to capture 
the main features of this activity which are relevant from 
a kinetic point of view. For hidden targets, the stim- 
uli emitted are very weak and the predator is very likely 
not to detect them. Then it may have to scan many 
times the same location before finding the target: we 
model this mechanism by a diffusion process. Note that 
in this description, the successive positions of the diffu- 
sive trajectory are not necessarily the very positions of 
the animal, but the focus points of the sense involved. 
This idea of a diffusion process (or more generally of 
random walk process) for sensory perception has already 
been suggested, in particular for vision[lj|, tactile sense 
or olfaction^. As focusing and processing the informa- 
tion received by sensory organs require a minimum time, 
the search phase can not be to short, which implies an 
upper bound /i max for the rate fa. Since the objective of 
phase 2 moves is to explore unscanned space, the most 
efficient solution, which is indeed observed generally, is 
a straight ballistic motion. Precisely, it has been ob- 
served that the turning angle after each pause is usually 
small [2! Hil Hii, Eof . Since the direction of successive bal- 
listic motions is strongly correlated for most of species 
(enjoying minimal memory skills), the small angles al- 
lows us to consider an effective 1-dimensional problem 
for both phases. 

We now evaluate the average time needed to find a 
target. Following the "closed cell approach" [f! Elf . our 
problem of an infinite space with uniform target density 



1/L can be reduced to the problem of a single target 
centered on a segment of size L with reflexive boundary 
conditions. Then, the instantaneous state of the searcher 
can be described by its position x on the segment and an 
index \x which specifies its motion: 

• /i = a (resp. (3) corresponds to a ballistic motion 
of velocity +v (resp. —v ) 

• /i = 7 (resp. 5) corresponds to a diffusive motion, 
switching only to a ballistic motion of velocity +v 
(resp. — v ) 

The mean first passage time at the target, starting from 
state (x,n) is denoted t(x,n). Using the Backward 
Chapman-Kolmogorov differential equation for the con- 
ditionnal densities [22I . we obtain the following system 
satisfied by the t(x, fx): 



dt(x, a) „ . . . 
v V ' + fa[t(x,y) 



D 



D 



t dtfx,f3) 

dx 
d t(x, 7) 

dx 
d 2 t{x,S) 

dx 2 



t(x, a)] 
fa[t{x,8)-t{x,(3)] 
fa [t(x,a) - t(x,j)] 
fa [t(x, (3) - t(x, 5)] 



= -1 



(1) 



Henceforth, we will consider the average search time 
(t) defined as the average of t(x, fj,) over the initial posi- 
tion x of the searcher, which is assumed to be uniformly 
distributed over the segment [—L,L], and over the na- 
ture of the initial motion, equally distributed over a and 
(3 (ballistic motions with velocities ±v). 

In the low density limit defined by L 3> 
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h ' V /: 

the system leads, after some calculation, to: 



(t) 
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rf! + 2 fa 
Vrfi + 4/i 



(2) 



where r = D/v 2 and 1/L is the target density. The linear 
dependence on L (the typical inter-target distance) en- 
sures that this combined strategy is much more efficient 
than a purely diffusive strategy^ which would scale like 
L 2 . The average search time (t) presents a single mini- 
mum with respect to fa and fa , defined by the following 
equations: 



fa = f: 



lmax 



/I + -/!/: 



1 ~ '2 j 1 " Z2/1 - 



(3) 



This minimum takes a simple form in two different 
regimes which lead to similar asymptotic: 

First, if /i max <ti 1/t, the optimal frequencies are such 
that fa = fa max and: 



fa = ( ^ ) ft 



1/3 



2/3 



(4) 
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In this regime, denoted S (for search), one has /i < fa: 
the predator spends more time searching than moving. 

Second, if /i max 3> 1/r the optimal frequencies are 
such that fx = f lmax and : 

/ 2 =(^y /3 /? /5 

In this regime, denoted by M. (for move), one has /i > fa 
and the predator spends more time moving. Note that 
the exponents 2/3 and 3/5 are numerically very close, and 
we do not expect to distinguish them experimentally. 

We now briefly comment on the structure of optimal 
trajectories. The threshold value 1/r has the meaning 
of an overlap limit: r is the duration for which the typ- 
ical distance covered is the same in both states 1 and 
2. In order to study the connexity of these trajecto- 
ries, we introduce the ratio R of the length scanned 
in phase 1 over the distance covered during a phase 
2 move R = y/D/fi x fa/ v. As for the S, one has 
R « (/it) 1 / 6 < 1 so that there is no overlap. It is 
noteworthy that in this regime, R can be small, which 
means that the overall scanned space is not connex but 
displays holes of unvisited space. Conversely, in regime 
JA, one has R « (/it) 1 ^ 10 > 1 and overlap may occur; 
nevertheless the power 1/10 ensures that R is never much 
larger than 1 and the optimal trajectory always explores 
unscanned spaces of significant size. 

We now compare our model with experimental data ex- 
tracted from O'Brien et al.0] and Kramer et al.Q which 
provide the average duration (and therefore its inverse, 
the rate) of search and motion phases for 18 different 
species (as various as planktivorous fish0], ground for- 
aging birds [TgL |24| or lizards[2^) performing a saltatory 
search behaviour. The corresponding rates f± range from 
0.1 to 100 Hz, with no systematic correlation with the 
animal size. Note that the prefactor involving r in equa- 
tion (@J and © depends a priori on the species and it 
seems difficult to evaluate it directly from observations, 
as D was defined in a phcnomcnological manner. Nev- 
ertheless, as we proceed to show, this characteristic time 
roughly assumes only two values. 

Using the exact relation J2J between fa, fa and r at 
the minimum, we compute the value of r for each species. 
The distribution of Log(r), given in figure[21 is quite un- 
expectedly bimodal. Remarkably, we observed that the 
first peak (around t=0.1s.) corresponds to foragers in 
the regime S and that the second one (around t=25s.) 
corresponds to foragers in regime M. Therefore the car- 
acteristic time t defined in our model appears as a tool 
characterizing distinctly two subclasses of foragers: a set 
S of animals in the regime S, and a set M in the regime 
M. . Since fluctuations of t are small within each set (the 
fluctuations of Log(r) are of order 1), we can perform a 
comparative analysis of these data within each set. 

The Log-Log plots of the fa and fa data of sets S and 
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FIG. 2: Log(r) distribution for saltatory search behaviours. 
The first peak (around t=0.1s.) corresponds to foragers in 
the regime S, the second one (around t=25s.) corresponds to 
foragers in regime M. 
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FIG. 3: Log-Log plot of experimental dataQ, of saltatory 
search behaviours, and their linear regression. 

M are shown in figure |3| Their linear regression shows 
that both sets are strongly correlated (with a coefficient 
r > 0.9) and that their slopes are in agreement with 
our theoretical prediction, namely 3/5 for M (experimen- 
tal slope: 0.7 ± 0.1) and 2/3 for S (experimental slope: 
0.6 ±0.1), taking account of the uncontrolled accuracy of 
experimental measurements. These results suggest that a 
wide variety of species indeed minimize their search time 
for preys according to a strategy described by our model, 
indicating that our diffusive modelling of the search phase 
is appropriate. This analysis also puts forward the hy- 
pothesis that the kinetics of the trajectories is a prevail- 
ing factor which is optimized by natural selection. 

Furthermore, using the extensive data of O'Brien et 
al.0 and Kramer et alQ, we have extended this anal- 
ysis to the case of animals performing a saltatory be- 
haviour observed in non specific search behaviours, such 
as fleeing, food carrying, or simple displacement. The 
distribution of Log(r) over these species is still bimodal 
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FIG. 4: Log(r) distribution for non specific saltatory be- 
haviours. 
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FIG. 5: Log-Log plot of experimental data0,0 of non specific 
saltatory behaviours and their linear regression. 

(see figure and highly similar to the case of figure 2, 
and leads to the definition of extended sets S' and M' as 
before. The corresponding Log-Log plots, presented in 
figure |5] show that the correlations between frequencies 
are maintained, and that the agreement with equations 
i jlj l and JSJ) is still significant: we find an experimental 
slope of 0.8 for both sets S' and M'. This suggests that 
all displacements of these species are conditioned by a 
search-like mechanism which could be either the search 
for a prey, or for a potential threat, as for example a hint 
of a nearby predator. Such an attitude could also be 
interpreted as a "behavioural economy" like argument, 
which would allow the animals to minimize the number 
of behaviours to learn. 

In summary, in this letter we have raised the ques- 
tion of determining the fastest strategy for finding a hid- 
den target. This question of general scope has been 
tackled through the example of animals searching for 
food, for which numerous experimental data are avail- 
able. We have proposed and solved analytically a two 
states stochastic model of target search, based on a diffu- 



sive modeling of the perception mechanism. This model 
provides a power law relationship between the character- 
istic times spent in each state, which satisfactorily fits 
experimental foraging data. Our findings suggest that 
in various behavioural contexts, saltatory animals adopt 
the intermittent motion which optimizes the search time 
of randomly hidden targets. We believe that this kind 
of search behaviour modeling could be extended to other 
situations, including human behaviours. 

Thanks to J-M Victor for fruitful discussions and to M. 
Barbi, R. Chitra, V. Fourcassie, M. Jardat, J-F Joanny , 
A. Lemarchand, E. Maisonhaute and P. Viot for critical 
reading of this manuscript. 
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